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Self-consistent simulations of the ultrafast electron dynamics in thin metal 
films are performed. A regime of nonlinear oscillations is observed, which 
corresponds to ballistic electrons bouncing back and forth against the film 
surfaces. When an oscillatory laser field is applied to the film, the field energy 
is partially absorbed by the electron gas. Maximum absorption occurs when 
the period of the external field matches the period of the nonlinear oscillations, 
which, for sodium films, lies in the infrared range. Possible experimental 
implementations are discussed. © 2008 Optical Society of America 

OCIS codes: 320.7110, 310.6860. 

The recent progress in the study of metallic nanostructures is mainly due to the development 
of ultrafast spectroscopy techniques, which allow the experimentalist to probe the electron 
dynamics on a femtosecond (and, more recently, attosecond) time scale. Typical "pump- 
probe" experiments involve perturbing the system via a first stronger pulse, followed by a 
second weaker pulse that acts as a diagnostic tool. By modulating the relative amplitude of 
the signals, as well as the delay between the pump and the probe, it is possible to assess 
with great precision the dynamical relaxation of the electron gas. 1-4 

In the present work, we focus on the ultrafast electron dynamics in thin metallic films. 
Several experiments have shown 1,2 that electron transport in thin metal films occurs on a 
femtosecond time scale and involves ballistic electrons traveling at the Fermi velocity of the 
metal Vp. These findings were corroborated by accurate numerical simulations, 5 which high- 
lighted a regime of slow nonlinear oscillations corresponding to ballistic electrons bouncing 
back and forth on the film surfaces. These oscillations were recently measured in transient re- 
flection experiments on thin gold films. 6 The existence of this regime prompted us to analyze 
the possibility of boosting energy absorption in the film by optically exciting the electron 
gas at the frequency of the nonlinear oscillations. 

In the rest of this Letter, time is normalized in units of the inverse plasmon frequency 
cjpg 1 , velocity in units of the Fermi speed Vf, and length in units of Lp = vp/uj pe . For alkali 
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metals we have L F = 0.59 (r s /a ) 1/2 A, u^ 1 = 1.33 x 1(T 2 (r s /a ) 3/2 fs, E F = 50.11 (r s /a )" 2 
eV and T F = 5.82 x 10 5 (r s /a ) K, where r s is the Wi gner-Seitz radius. We concentrate 
primarily on sodium films, for which r s = 4ao (oq = 0.529 A is the Bohr radius). 

We consider a system of electrons interacting via a Coulomb potential and confined within 
a slab of thickness L. The ion background is represented by a fixed density with soft edges, 
rii(x) = Hi [1 + exp [i\x\ — L/2)/o~i)] , where rTj = 3/(47rr^) is the ion density of the bulk 
metal and Oi -C L a diffuseness parameter. 7 In this jellium model, the self-consistent electro- 
static potential depends only on the coordinate normal to the surface (here noted x). Thus, 
the motion of an electron parallel to the surface of the film is completely decoupled from the 
motion normal to the surface and a one-dimensional (ID) model can be adopted. 

The electrons are initially prepared in a Fermi-Dirac equilibrium at finite (but small) 
temperature. They are subsequently excited by imposing a constant velocity shift Av x = 
0.08v F to the initial distribution. 7 This scenario is appropriate when no linear momentum is 
transferred parallel to the plane of the surface (i.e., q\\ = 0) and is relevant to the excitation 
of the film with optical pulses. 8 For q» = 0, only longitudinal modes (volume plasmon with 
u = uj pe ) can be excited. 

After the excitation is applied, the electron distribution function f e (x, v x , t) starts evolving 
in time according to the semiclassical Vlasov equation 

— +v— + — — — = (1) 

dt x dx m e dx dv x 

where m e is the electron mass and e denotes the absolute electron charge. The electrostatic 
potential is obtained self-consistently, at each instant, from Poisson's equation 

^ = -^-[n e (x,t) -m(x)] , (2) 

with n e = J f e dv x . As a reference case, we studied a sodium film with initial temperature 
T e = 0.0087V ~ 300 K, diffuseness parameter G{ = 0.3Lp, and thickness L = 50Lp — 59 A. 8 

The time evolution of the thermal E t h and center-of-mass E cm energies was analyzed 5 (Fig. 
[Q. During an initial rapidly-oscillating phase, E cm is almost entirely converted into thermal 
energy (Landau damping). After saturation, a slowly oscillating regime appears, with period 
equal to 50^;"^ pa 5.3 fs. This period is close to the time of flight of electrons traveling at the 
Fermi velocity and bouncing back and forth on the film surfaces (further details are provided 
in our previous work 5 ). 

The above nonlinear oscillations appear for all reasonable values of the physical parameters. 
Preliminary studies suggest that electron-electron collisions do not destroy this regime either, 
at least for relatively low excitation energies and short times. It is tempting, therefore, to 
investigate whether some kind of resonant absorption can be achieved when the system is 
externally excited at the same frequency of the nonlinear oscillations. 
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A similar scenario was investigated by Taguchi at al. 9 (building on an idea due to Brunei 10 ) 
in order to simulate the interaction of an argon cluster with a strong laser field (~ 10 15 — 
10 16 W/cm 2 ). In their simulations, the neutral cluster is quickly ionized by the laser field, 
which heats the electrons up to 10 eV. At these temperatures, the electrons behave classically 
and are initially described by a Maxwell-Boltzmann distribution. In that case, the electron 
transit velocity through the cluster is not clearly defined and depends on the intensity of the 
laser (indeed, in the classical case, there is no "natural" oscillatory regime like the one seen 
in Fig. 1). For a degenerate electron gas, the transit velocity is unambiguously given by the 
Fermi velocity of the metal and thus we expect an even neater resonance to occur. 

Our conjecture can be tested in the following way. At time uj pe t = 1000 (~ 106 fs for a 
sodium film), when the oscillatory regime is well established, we switch on a small external 
electric field, uniform in space and sinusoidal in time with period T: E cxt = E sin(27rt/T), 
where E is the (constant) field amplitude. The simulation is then continued with the external 
field on for another 4000c<j~ e 1 « 425 fs. This situation corresponds to a laser pulse that is 
switched on very quickly and lasts for a duration longer than 425 fs. 

For an electron transit velocity exactly equal to vp, we would expect resonance for a laser 
period T = 2L/vp (= 100 in units of uj^). The factor 2 comes from the fact that the electric 
field must keep the same sign during a transit from one surface to the other, and reverse sign 
during the "return" transit. We note that the resonance is expected to fall in the infrared 
(IR) domain. Indeed, for a laser period T = lOOo;^ 1 = 10.6 fs, the corresponding wavelength 
is A = 3.2 /j,m. 

The amplitude E of the laser field can be estimated by noting that the total energy of 
the laser pulse is U = (^EoEq) ctS, where c is the speed of light in vacuum, r is the pulse 
duration, and S is the surface of the laser spot. Typical values for IR lasers 11 are S = 0.01mm 2 
and U = l/i J, and by taking a pulse duration r = 400fs (similar to the duration used in the 
simulations), we obtain an electric field E = 4.3 x 10 8 V/m. In the numerical simulations, 
the electric field is normalized to E = m e VFU pe /e = 1.70 x 10 12 (r s /ao)~ 5 ^ 2 V/m, yielding 
E = 5.31 x 10 10 V/m for sodium films. Therefore, by taking a field amplitude Eq = 0.0LE, we 
get a dimensional value that is realistic for an IR laser pulse. This external field amplitude 
is an order of magnitude smaller than the self-consistent electric field present at the film 
surfaces. 

The results for the reference case (L = 50Lp) are shown in Fig. |2l where the electron 
thermal energy is plotted against time. We observe that the absorption is clearly enhanced 
for uj pe T = 106 and u pe T = 150, whereas for larger or smaller values virtually no energy is 
absorbed. We also verified that the resonance does not depend on the phase of the external 
oscillating field. 

The resonant period is close, but not exactly equal, to the predicted value uJ pe T = 100 
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and the resonance displays a certain broadness. The latter can be explained by noting that 
a certain dispersion exists in the electron velocities around vp, which generates a dispersion 
in the resonant period. If the period is T = 2L/v, then the resonance broadness should be 
\5T\ = (2L/v 2 ) 8v. In order to estimate the broadness, we plot, in Fig. El the variation of the 
velocity distribution at the center of the film: Sf(v x ) = \ f e (x = 0,v x ,LO pe t = 1000) — f e (x = 
0,v x ,t = 0)|. The distribution is indeed modified around the Fermi velocity, as expected (see 
also Fig. 3 in Ref. 5). Note that 5f is not symmetric around v x = 0, because the initial 
excitation was not symmetric either. On closer inspection, the peaks occur at a velocity 
slightly smaller (in absolute value) than vp, roughly \v x \ ~ 0.9vp. Their broadness can be 
estimated by assuming that a deviation of 2% (relative to the maximum f e — 1) is significant. 
Then, v varies in the interval 0.7t^ < v < I.Ivf (and equivalently for negative velocities), so 
that 5v ~ OAvp. This yields a broadness ui pe 5T ~ 50, with 90 < uJ pe T < 140. This estimate is 
compatible with the simulation results of Fig. 2, where the resonance has clearly disappeared 
at u pe T = 73 and 230. 

In order to test the robustness of this nonlinear resonance effect, we repeated the same 
numerical experiment with a thicker film, L = lOOLp (the initial evolution for this case is 
shown in our previous work 5 ). The resonant period is expected to scale linearly with the film 
thickness, and indeed we observed enhanced absorption for u pe T = 212 and 250 (Fig. 0]). 

In contrast, we observed that the resonance virtually disappears for thicker films, L = 
200Lp or larger. We interpret this result by noticing that the existence of the resonance 
depends on nonequilibrium electrons traveling coherently through the film. The phase space 
portraits of the electron distribution function (see Fig. 3 in our previous work 5 ) show a 
complex structure of traveling vortices. It is probable that, for thicknesses larger than a 
certain threshold, the necessary coherence is lost, so that the resonance cannot manifest 
itself. 

The resonance also disappears for very small amplitudes of the external field. For the 
reference case L = 50Lp, the resonance is still observed for Eq/E = 0.005, but no longer 
for Eq/E = 0.001. This may be related to the fact that the absorbed energy at resonance 
£/abs scales quadratically with the field amplitude: [/ a b s oc e 2 E 2 T 2 / m e (this formula becomes 
exact for the harmonic oscillator). For small fields, the resonance is thus very weak and other 
factors (e.g., Landau damping) can easily erase it completely. 

In summary, we have shown the existence of a nonlinear absorption regime in the electron 
dynamics of thin metal films. This effect is generic and should not, in principle, depend on 
the nature of the metal. The resonance occurs in the IR domain and should be accessible via 
experiments employing ultrafast laser sources with standard specifications. This absorption 
mechanism could be used as an optical diagnostic technique to determine, for instance, the 
thickness of the film, or to obtain information on the electronic distribution. 
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Fig. 2. 
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Fig. 4. 
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FIGURE CAPTIONS 



Fig. 1: Initial time evolution of the thermal and center-of-mass energies for a film thickness 
L = 50L F . 

Fig. 2: Time evolution of the thermal energy in the presence of an external electric 
field. The external field is switched on at u pe t = 1000. (a) u pe T = 27; (b) u pe T = 73; (c) 
ui pe T = 106; (d) u) pe T = 150; (e) u) pe T = 230. The results are for a film of thickness L = 50L F . 

Fig. 3: Variation of the electron velocity distribution with respect to the initial Fermi- 
Dirac equilibrium, at the center of the film, at time ui pe t = 1000. 

Fig. 4: Same as Fig. Elfor a film of thickness L = lOOLp- (a) u pe T = 90; (b) u pe T = 150; 
(c) u pe T = 212; (d) u pe T = 250; (e) u pe T = 290; (f) u pe T = 350. 



10 



